First detection of doubly deuterated hydrogen sulfide 
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ABSTRACT 

This work was carried out with using the Caltech Submillimeter Observatory 
and presents the observational study of HDS and D 2 S towards a sample of Class 



(N 
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sources, and dense cores. We report the first detection of doubly deuterated 
hydrogen sulfide (D2S) in two dense cores and analyze the chemistry of these 
molecules aiming to help understand the deuteration processes in the interstellar 
medium. The observed values of the D 2 S/HDS ratio, and upper limits, require 
an atomic D/H ratio in the accreting gas of 0.1 - 1. The study presented in 
this Letter supports the hypothesis that formaldehyde, methanol and hydrogen 
sulfide are formed on the grain surfaces, during the cold pre-stellar core phase, 



where the CO depleted gas has large atomic D/H ratios. The high values for 
the D/H ratios are consistent with the predictions of a recent gas-phase chemical 
model that includes Hjj" and its deuterated isotopomers, H2D + , D2H + and D3 



(Roberts et al. 2003). 
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1. Introduction 

Deuterium-bearing molecules have become the target of many observations in recent 
years and several models have been developed to account for them (e.g. Tielens 1983; Roberts 
& Millar 2000a,b). Twenty-six such molecules have been detected to date in interstellar 
clouds. The interest in this topic lies in the unusual chemistry at work in the cold regions 
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where CO is apparently strongly depleted (Bacmann et al. 2002, 2003) and in the general 
question of how much deuterium is trapped in the cold, dense phase of the ISM (Phillips & 
Vastel 2002). 

The study of doubly deuterated interstellar molecules has been very active in the last few 
years, since the surprising discovery of a large amount (~ 5%) of doubly deuterated formalde- 
hyde (D 2 CO) in the low mass protostar IRAS 16293-2422 (Ceccarelli et al. 1998; Loinard 
ct al. 2000). This is more than one order of magnitude higher than in Orion KL, where 
D 2 CO was first detected by Turner (1990). This first discovery was followed by many other 
studies which confirmed the presence of large amounts of D 2 CO (e.g. Ceccarelli et al. 2002) 
as well as doubly deuterated ammonia (ND 2 H) (e.g. Roueff et al. 2000; Loinard et al. 2001). 
Gas phase chemical models account relatively well for the observations of the abundances of 
singly deuterated molecules, but may not be able to completely reproduce the large deuter- 
ations observed for multiply deuterated molecules. The large deuterations could also be a 
product of active chemistry on the grain surfaces as predicted by Tielens (1983): deuteration 
during the mantle formation phase followed by evaporation of the mantles ices resulting from 
the heating of the newly formed star, with injection into the gas phase of the deuterated 
species (e.g. Ceccarelli et al. 2001). 

Recently, the detection of a triply deuterated molecule (ND 3 ) extended the observational 
framework. Until now, the possibility for detecting triply deuterated molecules seemed so 
remote that their lines were omitted in the spectroscopic catalogs for astrophysics. The 
ground-state rotational transition at 309.91 GHz of ND 3 has been detected with the CSO 
towards the Barnard 1 cloud (Lis et al. 2002) and the NGC 1333 IRAS 4A region (Van der 
Tak et al. 2002), with abundance ratios [ND 3 ]/[NH 3 ] ~ 10~ 3 and [ND 3 /H 2 ] ~ 10" 11 . The 
process at work, chemical fractionation, takes place in cold gas where the deuterium sub- 
stituted molecule is favored due to the lower vibrational uncertainty energy of the heavier 
molecule. The chemical process requires the presence of H3 which is usually limited by reac- 
tions with molecules such as CO. However, in these cold (~ 10 K) regions of the interstellar 
medium, molecules such as CO accrete onto grain surfaces to form mantles, so allowing the 
gas phase process of deuterium fractionation to take place (Bacmann et al. 2003). In time, 
the material of the mantle will suffer its own low temperature deuterium fractionation (Tie- 
lens 1983; Rodgers & Charnley 2002). Thus the grains may be important, both in enabling 
the gas phase process through depletion of reactive molecules and also by their own surface 
chemistry. For the latter to be effective, the deuterated molecules on the grain surface must 
be released, by some external energetic process, back into the gas phase, when they can be 
detected. The gas phase mechanism has recently received a boost from the recognition that 
inclusion of the more highly deuterated species D 2 H + and (Phillips & Vastel 2002) in the 
chemical calculations, increases at least by a factor of 2 the gas phase D/H ratio (Roberts 
et al. 2003). 



- 3- 



The identified most abundant components of the molecular mantles found on interstellar 
dust grains are H 2 and CO (for references, see Teixeira & Emerson 1999). The large range 
of condensation temperatures (see Nakagawa 1980) spanned by H 2 (T c ~ 90 K) and CO 
(T c ~ 16 K) suggests that molecules with intermediate condensation temperatures could 
also be present in the grain mantles. For example, H2S has a condensation temperature 
of ~ 39 K (Nakagawa 1980). This molecule and other sulfur-bearing molecules have been 
observed in the gas phase in cold, dark clouds where high-temperature processes should not 
be important (Minh et al. 1989; Irvine et al. 1989; Fuente et al. 1990): x(H 2 S) = 3 x 10~ 9 
at the SO peak of L134N, < 5 x 10~ 10 at the cyanopolyyne peak towards TMC1 and ~ 
7 x 10~ 10 towards the NH 3 peak position of TMC1. These values can be compared to the 
range 10 -11 - 10~ 10 obtained for steady-state gas-phase models (Millar & Herbst 1990; Millar 
et al. 1997). Therefore, it poses a problem for gas phase chemistry. The gas-phase reactions 
proposed for the formation of H 2 S are endothermic (see Section 3) and other mechanisms 
such as shock chemistry or high temperature chemistry are unlikely to be important in dark 
clouds. Grain surface reactions have then been suggested as an alternative pathway to H 2 S 
molecule formation (Duley et al. 1980). 

Because sulfur is chemically similar to oxygen, it is tempting to study this molecule to gain 
insight to the H 2 deuteration processes. Like H 2 0, it is possible that H 2 S originates on 
grain surfaces, where it was formed by accretion and hydrogenation of sulfur atoms during 
core formation (Tielens & Allamandola 1987). 

2. Observations and results 

The comparison between H 2 S, HDS and D 2 S would be useful for a complete study of the 
deuteration processes. In this letter, we report the observations of the HDS and D 2 S ground 
transitions, accessible at the Caltech Submillimeter Observatory (CSO). This observatory is 
a 10.4 meter single-dish antenna located on top of Mauna Kea in Hawaii. At 240 GHz, the 
telescope has an FWHM beam size of 31". The observations of these molecules were carried 
out between March 2002 and September 2002. Weather conditions were average with zenith 
opacities between 0.05 and 0.2 at 225 GHz. 

The frequency of the observed lines is 244555.58 MHz for HDS (l ,i - , ) with E u — 11.7 K, 
and 237903.6752 MHz for the ortho line of D 2 S (li,i - ,o) with E u = 11.4 K. These are 
the lowest lying transitions for the two molecules. The sample of sources observed includes 
three class (NGC1333-IRAS2; NGC1333-IRAS4A; IRAS16293-2422) and three dense cores 
(Barnard 1; the peak of emission in DCO + 3 — > 2 in the star-forming molecular cloud L1689N, 
Lis et al., 2002; the peak of emission in DCO + 3^2 near the young stellar object IRAS4A 
in the NGC1333 complex, Lis et al. in preparation). 
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The observation towards IRAS 16293-2422 was made in position switching mode with the 
reference position located at Sa = -180" and AS = 0" from the nominal center whose co- 
ordinates are quoted in Table 1. Positions 180" or 240" were observed using the chopping 
secondary mirror for the other sources. The main beam efficiency at 240 GHz is equal to 
70%. The spectra obtained towards these sources are presented in Figure 1 and the measured 
parameters in Table 2. 

The upper limits on the lines, when not detected, are derived following the relation: 



„ ^,^^/ r ^N 2 x binsize (km s x ) . , 

3 a (K = 3 x RMS (K x J ~ >- 1 

v ' v ' y Av (km s- 1 ) v ' 

and are reported in Table 2. From the observed line strengths given in Table 2, we can 
estimate the ratios between the D 2 S and HDS column densities in the observed sources. 
We estimated the optical depth from the observed line intensity using the formula: T m t> = 
\J v {T ex ) - J v (T bg )](l - e~ T ) where J V (T) = (hv /k) / (e hv ' kT - 1) is the radiation temperature 
of a blackbody at a temperature T, and T bg is the cosmic background temperature of 2.7 K. 
We found that the lines are optically thin, even at a low temperature of 10 K. The formula 
which connects the integrated intensity of a rotational transition with the number of emitting 
molecules is, for an optically thin transition and neglecting the background radiation: 



and, assuming the system is in LTE: 

The Einstein A coefficient is 1.3 x 10~ 5 s^ 1 for HDS (l ,i - ,o) (JPL catalog) and we 
computed a value of 5.0 x 10~ 5 s _1 for D 2 S (l^i - ,o)- Q(T ex ) is the partition function 
J2i9ie Ei/kTex equal to 3.1 for HDS (l ,i - , ), and 21.6 for D 2 S (1 M - 0)0 ) using an ex- 
citation temperature of 10 K. This temperature, taken as a kinetic temperature is that at 
which chemical fractionation should be efficient. The observed transitions are obviously not 
enough to calculate the excitation temperature in each source. Nevertheless, it is possible to 
estimate a reliable value for the excitation temperature using the kinetic temperature of the 
gas found in the literature. Menten et al. (1987) (respectively Bachiller et al. 1990) derived 
a kinetic temperature in L1689N (respectively Bl) of ~ 12 K based on their ammonia obser- 
vations. Blake et al. (1995) estimated T fc between 20 and 40 K for the core of IRAS4A based 
on their CS and H 2 CO observations. Ward-Thompson et al. (1996) found that T fc lies in the 
range of 15-20 K for IRAS2 from their infall model. Van Dishoeck et al. (1995) estimated a 
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rotational temperature of ~ 12 K from two HDS lines towards IRAS16293. In the three last 
sources, the Class sources, temperature gradients are known to be present, and the quoted 
estimates refer to the outer regions of the envelopes of these sources. Very likely the present 
HDS and D 2 S observations probe in fact those external regions too. In the following, we will 
assume that T ex is the same for HDS and D2S, namely that the 2 molecules originate in the 
same region. Luckily enough, the calculation of the N(D 2 S)/N(HDS) ratio does not depend 
so much on the excitation temperature in the 10 - 50 K range, as the energy levels compared 
are similar. 

D 2 S has para and ortho forms. The observed li 5 i - ,o transition is the lowest line in the 
ortho ladder. To derive the D 2 S/HDS ratio, the D 2 S column density was calculated for the 
ortho levels only, and then extrapolated to a total column density assuming an ortho to para 
ratio of 2. The derived D 2 S/HDS ratios have been computed using equations (1 - 3) and are 
reported in Table 3. 



3. Chemistry of D 2 S 

Hydrogen sulfide is traditionally believed to form on the grain surfaces, since gas phase 
reactions in cold gas are not efficient enough to account for the observed H 2 S abundances 
(e.g. Tielens & Allamandola 1987; Minh et al. 1989). In dense clouds sulfur is believed to 
be mostly neutral and initiates the chain reaction: 

S ^ HS+ H 3 S + ^ H 2 S (4) 

H 3 S + is thought to be the precursor to the H 2 S molecule, via a dissociative recom- 
bination. The reaction is important because all the normal SH n + + H 2 — > SH„ + i + + H 
reactions are highly endothermic (Millar et al. 1986) and cannot occur in cold interstellar 
clouds. Nevertheless, even at low temperatures, the rate coefficient of H 3 S + production is 
very low (~ 4 x 10~ 15 cm 3 s _1 at 10 K, Herbst et al. 1989), limiting the production of H 2 S 
molecules in cold clouds. On the contrary, sulfur is expected to easily react with H on grain 
surfaces and thereby form the volatile hydride H 2 S (Tielens & Allamandola 1987). In the 
same vein, formaldehyde and methanol are believed to be formed via active grain chemistry 
with successive hydrogenations of CO, stored on the grain mantles, and eventually released 
into the gas, for instance during the collapse phase, when the heating of newly formed pro- 
tostars evaporates the CO-rich ices. 

If the deuteration of the three ices (H 2 S, H 2 CO and CH 3 OH) occurs when the ices are 
formed and if it is dominated by pure statistics (i.e. by the random distribution of D atoms), 
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the deuteration ratios would assume the following values (detailed discussion in Rodgers & 
Charnley 2002): F(H 2 S) = (HDS/H 2 S) 2 /(D 2 S/H 2 S) = 4, F(H 2 CO) = (HDCO/H 2 CO) 2 /(D 2 CO/H 2 CO) 
= 4, F(CH 3 OH) = (CH 2 DOH/CH 3 OH) 2 /(CHD 2 OH/CH 3 OH) = 3. Observed values differ- 
ent from these would indicate either different chemical pathways or chemical reactions on 
the grain surfaces which favor D over H (or vice- versa) substitutions. Note that the modified 
rate approach and Monte Carlo method by Caselli et al. (2002) to explain deuteration frac- 
tionation on interstellar grains indeed confirm the above values -assuming that the chemistry 
of S is similar to that of O on the grain surfaces. 

Our sample includes IRAS16293-2422, where the singly and doubly deuterated forms of 
both formaldehyde and methanol have been observed (Ceccarelli et al. f 998; Loinard et al. 
2000; Parise et al. 2002). The observed values are: F(H 2 CO) = 0.4 and F(CH 3 OH) = 4. 
Combining the present D 2 S/HDS upper limit (Table 3) with the HDS/H 2 S value found by 
Van Dishoeck et al. (1995, HDS/H 2 S~ 0.1) gives: F(H 2 S) > 0.5. It seems hence that while 
methanol and hydrogen sulfide are consistent with the hypothesis of the random distribution 
of D's and formation on the grain surfaces, formaldehyde is definitively off. The use of the F 
factors has the great advantage to highlighting what chemical mechanism is responsible for 
the observed deuteration (Rodgers & Charnley 2002), but obviously overlooks the details of 
it. 

We computed in Table 3 the column densities for the HDS molecule. Assuming that the 
HDS/H 2 S is around 10%, the value found for formaldehyde and other molecules presenting 
a large degree of deuteration, one would obtain H 2 S abundances in the range (0.2 - 2) x 
10~ 9 . This value has been computed with an H 2 column density about 2.8 x 10 23 cm -2 
for L1689N(DCO+) (Lis et al. 2002), 2 x 10 23 cm" 2 for IRAS16293 (Van Dishoeck et al. 
1995), 1.3 x 10 23 cm" 2 for Bl (Hirano et al. 1999) and 3.1 x 10 23 cm" 2 for IRAS4A (Van 
der Tak et al. 2002). We don't know yet the actual value for the H 2 column density towards 
IRAS4A(DCO + ). We can say that the computed H 2 S abundance seems higher than the 
range of 10~ n - 10~ 10 estimated for steady-state gas-phase models (Millar & Herbst 1990; 
Millar et al. 1997), and it would be a strong evidence that grains, and not gas, are producing 
H 2 S. Of course, this is a rough estimate, and observations for the H 2 S molecule have already 
been planned in order to derive more accurate values. 

The diagonal line in Fig. 2 shows the predicted D 2 S/HDS as function of the atomic D/H 
ratio in the accreting gas, assuming that D 2 S/HDS is equal to D 2 0/HDO, as computed by 
Caselli et al. (2002). The argument that the surface chemistries of sulfur and oxygen are 
similar can be challenged by the fact that D 2 has not yet been detected in the interstellar 
medium. But the search for D 2 has not yet been carried out towards the same objects 
where D 2 S has been detected, and it will be interesting indeed to check whether the two 
molecules, D 2 S and D 2 0, have very different abundances with respect to their respective 
main isotopomers. In the model presented by Caselli et al. (2002), surface chemistry can 



-7- 



convert a high atomic D/H ratio in the gas to high abundance ratios for singly, doubly or 
triply deuterated isotopomers, especially at high gas densities (10 5 cm -3 ) and low temper- 
atures (10 K). The observed values in Bl and NGC1333-IRAS4A(DCO+) and upper limits 
of Table 3 are consistent with previous determinations of the required atomic D/H ratio in 
the accreting gas of ~ 0.2 - 0.3 found in the case of IRAS16293-2422 (Parise et al. 2002). 
We conclude that, based on the observations so far available, hydrogen sulfide, formaldehyde 
and methanol are very likely formed on the grain surfaces, and the deuteration occurs during 
the formation of the mantles. What remained problematic to explain up to very recently 
is the relatively large value of the atomic D/H ratio required to achieve the observed large 
deuteration fractionation in the three molecules. In fact until very recently, no gas-phase 
chemical model predicted such large values of atomic D/H (e.g. Roberts & Millar 2000a,b). 
However, the current year has seen the detection of an extremely large abundance of H 2 D + in 
the pre-stellar core L1544 (Caselli et al. 2003), 10 times larger than prediction from previous 
chemical model. This detection, together with the previous observations of the doubly and 
triply deuterated molecules, has stimulated a radical revision of the chemical gas phase mod- 
els. Phillips & Vastel (2002) pointed out that the Hg chemistry should include D 2 H + and 
D3 , as well as the expected H 2 D + . Roberts et al. (2003) showed that inclusion of the doubly 
and triply deuterated forms of H3 would increase the atomic D/H ratio in CO depleted gas 
to the levels 0.1 - 1, required by the D 2 CO, CHD 2 OH and the present D 2 S observations. In 
conclusion, this work supports the hypothesis that formaldehyde, methanol and hydrogen 
sulfide are indeed formed on the grain surfaces, during the cold pre-stellar core phase, where 
the CO depleted gas possesses large atomic D/H ratios. 
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Fig. 1. — Spectra acquired with CSO of transitions of HDS and D2S towards a sample of 
class and dense cores. X-axis represents the velocity in the local standard of rest and 
Y-axis represents the main beam temperature in Kelvins. 
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Fig. 2. — D 2 S/HDS abundance ratio plotted as a function of the D/H abundance (cm -3 ), for 
high density conditions at 10 K using a modified rate approach of formation on interstellar 
grain surfaces (see Caselli et al. 2002). Bl and NGC1333-IRAS4A are plotted at the same 
D 2 S/HDS ratio, but Bl is a detection and NGC1333-IRAS4A an upper limit only. Upper 
limits are represented with arrows. 
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Table 1. Source coordinates. 



Source 


RA (2000) 


DEC (2000) 


NGC1333-IRAS2 


03h28m55.4s 


31°14'35.1" 


NGC1333-IRAS4A (DCO+) a 


03h29ml2.1s 


31°13"26.4" 


IRAS16293 


16h32m22.6s 


-24°28"33.0 // 


Bl 


03h33s20.8s 


31°07'34.4" 


NGC1333-IRAS4A 


03h29ml0.3s 


31°13'32.3" 


L1689N (DCO+) b 


16h32m28.6s 


-24°29'2.7" 



a This source corresponds to a position 23" East and 6" 
South of the IRAS 4A position, and is where the DCO + 3 — > 
2 emission peaks (Lis et al., in preparation) 

b This position corresponds to the DCO + emission peak 
observed by Lis et al. (2002) in the star forming molecular 
cloud L1689N 
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Table 2. Line parameters. 





HDS 


(lo,i-0 o ,o) 


D 2 S 


(li,i-M 




(244 555.58 MHz) 


(237 903.6752 MHz) 


Object 




Av 




Av 




(K) 


(km/s) 


(K) 


(km/s) 


NGC1333-IRAS2 


< 0.05 




<0.03 




NGC1333-IRAS4A (DCO+) 


0.31 


1.4 


0.08 


1.8 


IRAS16293 


0.13 


1.4 


< 0.1 




Bl 


0.22 


1.1 


0.07 


0.6 


NGC1333-IRAS4A 


0.17 


1.7 


< 0.05 




L1689N (DCO+) 


0.15 


1.5 


< 0.06 
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Table 3. D 2 S/HDS ratio and HDS column densities for an excitation temperature 

between 10 K and 50 K. 



N(D 2 S)/N(HDS) N(HDS) 



Object 


(%) 


(10 13 


cm" 2 ) 


NGC1333-IRAS4A (DCO+) 


12.0 


1.2 


- 4.8 


IRAS16293 


< 26.4 


0.5 


- 2.0 


Bl 


10.8 


0.7 


- 2.6 


NGC1333-IRAS4A 


< 10.8 


0.8 


- 3.1 


L1689N (DCO+) 


< 14.4 


0.6 


- 2.5 



